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Within the Schwinger-Keldysh formalism we derive a Ginzburg-Landau theory for the Bose- 
Hubbard model which describes the real-time dynamics of the complex order parameter field. Ana- 
lyzing the excitations in the vicinity of the quantum phase transition it turns out that particle/hole 
dispersions in the Mott phase map continuously onto corresponding amplitude/phase excitations in 
the superfluid phase, which have been detected recently by Bragg spectroscopy measurements. 



I. INTRODUCTION 



Within the last decade ultracold atoms in optical lat- 
tices [U, [H have become a standard tool for studying 
quantum-statistical many-body effects. Due to their high 
tunability, these systems represent an almost perfect test 
ground for a large variety of solid-state models. In partic- 
ular, the experimental observation of the seminal quan- 
tum phase transition from the Mott insulating (MI) to 
the superfluid (SF) phase, exhibited by a single-band 
Bose-Hubbard (BH) system of spinless or spin-polarized 
bosons, has recently attracted a lot of attention Al- 
though the occurrence of this quantum phase transition 
is evident from the momentum distributions of time-of- 
flight measurements, its precise location cannot be deter- 
mined from them. Recently, however, more detailed in- 
formation about the collective excitations of this system 
could also be achieved by exciting the system via lattice 
modulation [f| or by Bragg spectroscopy Deep in 

the SF phase, the observed gapless excitation spectrum 
can be well understood within a Bogoliubov theory 
Approaching the phase boundary, a time-dependent dy- 
namic Gutzwiller calculation [7], a slave-boson method 
[l(]| |. and a random-phase approximation [llj have pre- 
dicted an additional SF gapped mode, which recently 
could be confirmed experimentally in the strongly inter- 
acting regime 7]. An open question is the fate of this 
mode in the weakly interacting Bogoliubov limit. When 
the MI phase is reached, both SF modes turn contin- 
uously into particle and hole excitations which are also 
found by mean- field theory [§] . Due to their finite energy 
gaps, they characterize the insulating phase in a unique 
way. 

A field-theoretic ansatz describing the system in both 
the insulating and the superfluid regime has first been 
considered in Ref . [l2j , where an effective action has been 
obtained via two successive Hubbard-Stratonovich trans- 
formations. The same action can also be obtained by the 
Ginzburg-Landau approach developed in Refs. [H, [13], 
which is technically based on resumming a perturbative 
hopping expansion. As is explicitly shown in Refs. [l3l - 
[l7j . the practical advantage of this non-perturbative ap- 
proach is that it provides a generalization of mean-field 
theory by taking into account higher hopping orders in 



a straight-forward systematic way. However, the Lan- 
dau theory Ref. [13| is restricted to a static descrip- 
tion of the SF-MI transition at zero temperature. A 
finite-temperature Ginzburg-Landau theory developed in 
Ref. fl4| yields also dynamic results via an analytic con- 
tinuation, but is restricted to near-equilibrium situa- 
tions. Since time-resolved measurements jl8| have be- 
come possible within the last years, a real-time descrip- 
tion of quantum systems is desirable. Therefore, we fol- 
low Ref. [1911 and modify the imaginary-time approach 
from Ref. [lj] by converting it to real time, thus requiring 
techniques which were first introduced by L.V. Keldysh 
and J. Schwinger 

To this end we start with a brief introduction to the 
so-called Schwinger-Keldysh formalism of real-time per- 
turbation theory at finite temperature in Section UT] The 
actual Ginzburg-Landau theory for the underlying Bosc- 
Hubbard model is then derived by adapting the non- 
perturbative imaginary-time calculation of Ref. [l4[ to 
the Keldysh space. Thus, we combine a perturbative 
hopping expansion of the free energy in Section IIIII with 
a subsequent resummation due to a Legendre transfor- 
mation in Section IIVI and determine the resulting ef- 
fective action. Its extremization leads to equations of 
motion for the order parameter fields which are solved 
in Section [V] yielding a plethora of physical results like 
the phase boundary and the excitation spectra in Sec- 
tion IVII Furthermore, we compare in Section EH] the 
present Schwinger-Keldysh calculation with the corre- 
sponding imaginary-time approach of Ref. fl4j , revealing 
an unexpected mismatch between both descriptions. Fi- 
nally, a summary of the present work, which is restricted 
to equilibrium physics, and an outlook, how to generalize 
it to non-equilibrium situations, is given in Section [Villi 



II. SCHWINGER-KELDYSH FORMALISM 

In order to introduce the Schwinger-Keldysh formal- 
ism, we consider the Green functions of a many-body 
system. In general, these are averages over a product 
of creation operators a] and annihilation operators a,j, 
where the index i collects the degrees of freedom of the 
system. In our case, it just denotes the respective lattice 
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site. At zero temperature, the averaging process is purely 
quantum mechanical, whereas for finite temperature also 
the thermal fluctuations are taken into account. In equi- 
librium systems, i.e. those with a Hamiltonian constant 
in time, introducing imaginary-time variables puts the 
thermal averaging formally on an equal footing with the 
dynamic evolution of a time-dependent system at zero 
temperature. 

For both cases, perturbative treatments are well 
known, which arc based on reducing the non-trivial real 
or imaginary dynamics of the operators in the Heisen- 
berg picture to the trivial dynamics of operators in the 
Dirac picture, which is determined by some solvable part 
Ho of the full Hamiltonian H. The remaining part Hi 
then enters the time-evolution operator U(t,i ) which 
acts on the states. In the case of a real-time evolution, for 

instance, it reads V(t,t ) = Texp — | f f dtiHiD(fi) 
where T is the time-ordering operator. 

In the zero-temperature formalism (ZTF), the Green 
functions are defined as ground-state expectation values 
of time-ordered operator products in the Heisenberg pic- 



ture, i.e. (T 



notes the expectation value with respect to the ground- 
state of the Hamiltonian H. Applying the translation 
formula from the Heisenberg to the Dirac picture to this 
expression, i.e. an(t) = \J(to,t)&n(i)XJ(t,to), where to is 
the initial time, we find for each operator a forward and a 
backward time-evolution. But due to the time-ordering, 
these various pieces can be straightened to one single 
forward evolution from to to the largest time within the 
operator product and a subsequent backward evolution. 
The latter, however, cancels due to a connection between 
the ground-states of the full and the unperturbed sys- 
tem, which is provided by the Gell-Mann-Low theorem 
[22j . In this way, we do not only get rid of the backward 
evolution, but also reduce the expectation value to the 
one of the known, unperturbed ground state. Finally we 

U(oo, -oo)aj iD (ti) • • • a,i n T>(t n ) ) . , where the 

L J / H 

non-trivial part of the averaging process enters only via 
a forward evolution along a straight time path, which 
has been extended to infinity. In this way, the pertur- 
bation theory is reduced to a Taylor expansion of the 
time-evolution operator U(oo, —00). 

In the following, however, we also consider finite tem- 
peratures, so in the Green functions the average (») H now 
denotes a thermal average with respect to the equilibrium 

of H, i.e. Tr (• e"'™) /Tr (e^ 6 ) , where /3 = l/k B T is 

as usual the inverse temperature. If there is no evolution 
along the real-time axis, then the scheme of the ZTF 
described above can be taken over to an artificially in- 
troduced imaginary-time axis. Evolution along this axis 
somehow imitates the thermal averaging process, so the 
backward evolution in imaginary dynamics of the opera- 
tor product can directly be cancelled against the pertur- 
bative part of the thermal density matrix e~^. In this 
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way, the imaginary-time formalism (ITF) also puts all 
non-trivial Hi contributions of the Green functions into 
a straight time evolution. 

The situation becomes more complicated for a time- 
dependent many-body system at finite temperature 
where the backward evolution along the real-time axis 
cannot be canceled by the Gell-Mann-Low theorem as in 
the ZTF. To realize the consequences of this, let us con- 
sider for simplicity the product 2 n(t 2 )Om(ti) of two 
arbitrary operators in the Heisenberg picture. If t\ > 
t2, the translation from the Heisenberg into the Dirac 
picture yields U(t , t 2 )6 2D (t 2 )XJ(t 2 , ii)6i D (*i)U(*i, to), 
which represents a time-evolution along a closed contour 
as depicted in Fig. Q] a) with G^d^) being located on 
the backward path of the contour. For the opposite case 
t 2 > tx, the closed contour extends from to to t 2 , and the 
condition to maintain the order of the operators now is 
that Oid(^i) appears on the forward path. In both cases 
it plays obviously no role, which part of the contour is 
assigned to the last operator. 

By providing each operator with a corresponding path 
index P = ±. where + stands for the forward and 
— for the backward path, the positions of the opera- 
tors on the contour can be fixed. With this we are 
able to define a contour-ordering operator T c , which 
brings the +operators in time order, while acting as an 
anti-time ordering operator A on the —operators. Such 
contour-ordered operator products are the most natu- 
ral generalization of the time-ordered products consid- 
ered in the ZTF. Thus the relevant Green functions are 



given by quantities like ( T, 



But even after transforming into the Dirac picture, 



H 

i.e. 



where 



t > is the biggest time appearing in the operator prod- 
uct, the perturbative part of the Hamiltonian is not com- 
pletely isolated within the time-evolution operator, since 
it still appears in the thermal average. We might man- 
age this by extending the time-evolution contour to the 
complex plane as shown in Fig. [T]b), yielding a third 
path. 

However, it is widely believed in the literature that, if 
we perform the initial time limit to —00, the imaginary 
part of that contour can be neglected for initially uncor- 
rected systems [H, [24|. We will follow this tradition 
and discuss its limitations at the end of this paper. With 
this simplification and by extending the contour to the 
infinite future, we get the time-evolution contour from 
Fig.Q]c), encircling completely the real-time axis. It is a 
closed path, and we, therefore, refer to this formalism as 
the closed time path formalism (CTPF). 

Then the contour-ordered Green functions can be de- 
fined as 

G|f } } ({t}) = i n+m ~ '{t c [s-tS a£(*i) • • • af;(t n ) 

xaf-^+i) • • • <; + m mt (w)l ) . , (1) 

J / Ho 

where we have waived for simplicity the index for the 
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FIG. 1: Contours of time evolution: a) Schwinger contour, b) 
interaction contour, c) Keldysh contour. 



Dirac picture, as we will do so for the rest of this 
paper. With the curly brackets {X}, we abbreviated 
the union of variables X± , • • • , X n+m . In this expres- 
sion, the perturbation enters only in the time evolution 



S' S = T c exp 



Br 



a+(t)- ar w]}. 

The Green functions from Eq. ([1]) can systemati- 
cally be derived from a generating functional. To this 
end, we have to introduce artificial currents jf(t) and 
their complex conjugate jf(t), which linearly couple 
to the corresponding creation and annihilation opera- 
tors af^(i) and af (t). With such a source term, the 
new Hamiltonian reads H p [j, j] (t) = H p (i) + H p (t) + 

J2i if (*)&f (*) +jf (*)*f(*) • Assigning the source 
term to the perturbative part, we have to substitute S by 
S\j+,j+] = Texp (-U-o. M&tW+^i [3tmP(t)+ 
c.c] and similarly for the backward evolution opera- 
tor §t. The generating functional is then defined as 



Tc{St[7-,r]S[j- + ,j + ]}) , (2) 

Ho 



yielding the Green functions as its functional derivatives 

3=0 



■■•i n ,i n +i---i n + v 



S n+m Z[jJ] 



Sjh • • • Sji n Sj in+1 ■■■Sji n ^ 



(3) 



3=0 



where both path indices and temporal variables have 
been waived for simplicity. 

Path-ordered quantities, though naturally arising 
within the Keldysh formalism, lack a clear physical in- 
terpretation, so it is convenient to define the linear com- 
binations 



X E (t) = [X+(t) + X-(t)]/V2, 



X*{t) = [X + (t)-X-(t)]/V2 



(4) 
(5) 



A good discussion of different choices for a suitable basis 
can be found in Ref. (25|. To see how the path-ordered 
Green functions given in Eq. ([I} transform under this so- 
called Keldysh rotation, it is most instructive to consider 
the 2-point functions. Writing them in form of a 2x2 
matrix, the transformation is given by 



G± + (t;t>) G±-(t,t>) 



GfjW) 



(6) 



where G R ' A are the well known retarded/advanced Green 
function and A is the anticommutator function. In the 
case of bosons, they are defined as (2(| 



Gf 3 {t-t<) 
G^t-t') = - 
Aij(t;t') = i 



i6{t-t')( [a iH (i),at H (t') 
-i9(t-t')( [a iH (t),at H (t') 



(7) 
(8) 
(9) 



Here, [•, •] + denotes the anticommutator, while the com- 
mutator is given by [•,•], and 9 stands for the Heavi- 
side function. To see how these quantities arise from the 
Keldysh rotation in Eq. ^ , we must note that the defini- 
tions in Eqs. ©-© can be rewritten in terms of contour- 
ordered products T c [a^ I (t)aJ I /(t')] . It then turns out 
that these can be summarized to one single product of 
a E and a A operators or their hermitian conjugates, re- 
spectively. We find: 



Gf j (t ] t') = i(f c [kf s (t)kfi(t')] 



G^(*;f) = i<Tc[ag I (t)aS(0] 



A ij (t;t') = i{T c [kf H (t)k%(t')] 



(10) 

(11) 
(12) 



Since in the next section we are going to calculate the 
path-ordered Green functions perturbatively, we stress 
that these relations do not only hold for the exact 
Green functions, but also in any order of perturbation 
theory (27| . Furthermore, it can be shown for arbi- 
trary n-point functions that the Keldysh rotation always 
yields, amongst other functions, the retarded and ad- 
vanced Green function. For practical calculations it is 
also relevant to notice that operator products of the kind 
T c O A • • • A are identically zero [23]], so amongst the 
Keldysh-rotated n-point functions, we always find a van- 
ishing one. 



III. PERTURBATION THEORY 

In the previous section we have introduced the gen- 
eral formalism of separating the non-trivial real-time dy- 
namics and thermodynamics from the rest of a Green 
function by performing a translation from the Heisen- 
berg to the Dirac picture. In view of concrete applica- 
tions, we must now define the respective unperturbed 
and perturbed parts of the considered Hamiltonian. 

Bosons in an optical lattice are described by the BH- 
Hamiltonian [3, 0] given as Hbh = Ho + Hi , where we 
have introduced a decomposition into the local part 



u 



a^aj 



(at&i - l) 



aj 



and the non-local hopping term 



Hi 



Jij ^3 ■ 

i,3 



(13) 



(14) 
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Here a, (aj) denotes the bosonic annihilation (creation) 
operator at lattice site i, [i the chemical potential, U the 
on-site interaction parameter, and Jij the hopping matrix 
element being equal to J > for nearest neighbors only. 

The main feature of this Hamiltonian is a quantum 
phase transition from a MI phase to a SF phase due to 
the competition between kinetic energy, i.e. hopping be- 
tween sites, and local on-site interactions. Therefore, it is 
natural to consider the hopping as a perturbation, if we 
are interested in the MI phase, whereas in the SF phase 
the interaction represents the perturbation. However, a 
perturbative description of the critical behavior in be- 
tween both regimes would then be doomed to fail. But 
we are in the lucky situation that a resummed hopping 
expansion can also be considered as a 1 //^-expansion of a 
D-dimensional system and, thus, might allow for a proper 
description of the quantum phase transition. To this 
end we note that, in the presence of a condensate, the 
hopping parameter J must be rescaled according to the 
dimensional scaling law J — > J / D in order to have a fi- 
nite energy in infinite dimensions [28j . This means that 
within a hopping expansion all nth order hopping loops 
are suppressed by a factor 1/D n . By resumming the 1- 
particle irreducible contributions up to the nth hopping 
order, we therefore get an effective (l/D)-expansion up 
to the (n — l)th order. Although we restrict ourselves 
in this paper to the lowest order n = 1, we will get a 
theory which is exact for infinite dimensions or infinite- 
range hopping [291] . In this way we find for D > 2 that 
the quantum phase transition can be well described in a 
quantitative way. 

Finally we must add a source term to the Hamiltonian, 
yielding H[j, ]}(t) = H BH + H s (t), with 



fts[7,i](*) = E[^(*)M*) 



(15) 



A similar source term has already been introduced for 
defining a generating functional in Eq. @, but there it 
has only been a technical tool. Although we will set the 
currents to zero in the end, here this term acquires a 
physical meaning, since it explicitly breaks the underly- 
ing U(l) symmetry of the original BH-Hamiltonian. Only 
this symmetry-breaking makes it possible to have a non- 



vanishing condensate amplitude (a,i(t)), which character- 
izes the SF phase and can be taken as an order parameter 
field ^i(t) for describing the MI-SF transition within a 
Ginzburg-Landau theory. From the definition (1151) it is 
clear that the currents ji(t) and ji(t) are the conjugate 
variables of the order parameter fields ^i(t) = (aj(i)) 
and ^i(t) — (ai(i)), now appearing explicitly within the 
Hamiltonian H[j, j](i). 

Ginzburg-Landau theories are in general based on a 
fourth-order expansion of the thermodynamic potential 
in the order parameter field [30] . In order to derive it for 
the underlying Bose-Hubbard model, we will start from 
the free energy. For small order parameter fields we can 
assume that also the influence of the symmetry-breaking 
currents is small, so we can write down the free energy 
as a power series with respect to these currents ji (t) and 
ji(t). It is then possible to replace the currents via a 
Legendre transformation by the order parameter fields 
%(t) = (a|(£)) and <F;(i) = (a;(i)), so that we finally 
end up with the so-called effective action that serves as 
a Ginzburg-Landau functional. 

Noting that the generating functional Z[j + , j~, j + , j~] 
defined in Eq. ^ has a similar structure as the partition 
function in statistical mechanics, we get an analogue of 
the free energy by taking the logarithm: 



J r [j + ,j ,j + ,j~ 



iln Z\j+,j-,j + ,j-}. (16) 



As described in Section UH the functional Z is given by 
a forward and a backward evolution in time, where the 
operators and currents appearing on both time paths 
have to be distinguished from each other. Before start- 
ing a perturbative expansion of this functional, however, 
it turns out to be feasible to perform a Keldysh ro- 
tation, which mixes the quantities on the forward and 
backward path as defined by Eqs. (|4]) and ([5]). Un- 
der this rotation, the source term transforms as jf&l — 



have a, + ^at — a, ^a ■ 



j i a^ T = jf &f * + jj A af \ and for the hopping term we 
a^a? + af ^a A . It turns out to be 
helpful to introduce the vector quantities ji = (jf^,jf) T 
and aj = (af^,a^). Then the generating functional 
reads 



/ ■ />oo 

H J] = -i In ( Tc exp ( - - J dt{ J2 
where a 1 is the Pauli matrix I31J. 



Instead of straightforwardly expanding this functional in 
the hopping parameter J and in the vector currents j and 
j, we first note that the logarithm has the nice property 
to make .F[j,j] an extensive quantity. Thus the linked- 



,a] ( r 1 a l + ^[j l (i)-a l (t)+c.c.]})\ , (17) 

I 

cluster theorem [H, [33| applies: It states that, whereas 
the expansion of Z is built up by the Green functions, 
which have been defined in the path-ordered basis in 
Eq. (fTJ), the functional T can be expanded in a series 



5 



of connected Green functions or cumulants C: 



J = in Eq. JT]). Explicitly, they read: 



3=0 



.(18) 



3=0 



For simplicity, we have suppressed here again both the 
time variables and the Keldysh indices. 

These cumulants are much simpler objects than the 
Green functions: Whereas roughly speaking the Green 
functions may describe any particle creation and anni- 
hilation processes, the cumulants decompose them into 
their independent contributions. What is meant by this, 
can be illustrated by considering only the unperturbed 
system Ho in (fT5|) . Here, each lattice site i can be con- 
sidered as an independent system and we could define 
generating functionals -2rf°' ) [ji , ji] in each subsystem. We 
would find that the Green functions of the whole system 
could be derived from Z(°)[j,j] = JliZ- ■ [jijj]. Whereas 
an expansion of each .Z^ -* [ji, ji] in the currents would be 
completely local, the expansion of would con- 

tain highly non-local objects as well. Nevertheless it is 
clear, that these non-local Green functions arise as prod- 
ucts of local Green functions, and it must therefore be 
possible to find a decomposition. This mixing is cir- 
cumvented from the beginning, if we consider instead 

•?" (0) [jJ] = -lY^i^^i [jiJi]) where the logarithm has 
turned the product of subsystems into a sum of them. 

In the system with hopping, of course, the situation is 
not that simple, but if we expand not only in the 

currents, but also in the hopping parameter J, the expan- 
sion still contains only the local cumulants of the unper- 
turbed system, and the hopping simply appears as a link 
between two cumulants. For an expansion up to fourth 
order in the currents, we will need the unperturbed cu- 

mulants Cf lK2(0) (t 1 ;t 2 ) and Cf lK ^ m {h^MM), 
where K — £, A denotes the respective Keldysh index. 
All the cumulants have only one site index due to their lo- 
cality. Note that cumulants describing an unequal num- 
ber of annihilation and creation processes vanish in the 
unperturbed system, since Ho commutes with the local 
number operator fij = aja»j. Their relations to the un- 
perturbed Green functions can be derived from Eqs. ([3]) 
and l|18p and read explicitly 



a. 



KiKa(O) 



(h;t 2 ) = G 



KiK 2 (0) 



(h;t 2 ), 



and 



{^1,^2,^3,14,) 
_ ^KxKaKsK^O)/. s .^K 1 K 3 (0), s 



(19) 



(20) 



xCf 2K ^ 0) (i 2 ;< 4 )- i Cf lK ^ 0) (ti;t4)Cf 2K3(0) (i 2 ;f 3 ). 
The unperturbed Green functions are obtained by setting 







} 


Tr- 


' e -/3Ho| 



and 



r ,K 1 K 2 K 3 K 4 (0)/. , . i \ 
Tr- 



(21) 
(22) 



Tr 



{e-^Hoj 



As argued at the end of Section HH we have G AA = 

C AA = gAAAA = C AAAA = from which follows that 

c k iK2 k 3 k 4 = 6 tJ 6 jk 6 kl Gf/ k f 2K3K4 , if three of the four 
Keldysh indices Ki,K 2 ,K 3 ,K4 are A. 

It is feasible to interpret the 2n-point cumulants as 
tensors of rank 2n in a two-dimensional space, which ac- 
counts for the doubled time degree of freedom. We denote 
these tensors by C,(ti, • • ■ , t n ; t n +i, ■•• , i 2 «)- The most 
important objects within the Schwinger-Keldysh formal- 
ism are the two-point cumulants, which form a 2x2 ma- 
trix: 



Cj(ii, ; i 2 ) 



C^°\t u t 2 ) Gf A(0) (^,t 2 ) 







(23) 



A comparison with Eq. ^ together with Eq. (fTi?)) yields 
that the upper (lower) off-diagonal element are the re- 

R /A 

tarded (advanced) Green functions G u (ii;f 2 ). 

With the definition of vector operators and vector cur- 
rents as in Eq. (|17p and the definition of tensor cumu- 
lants, we can formally circumvent using explicitly path 
or Keldysh indices. Within a diagrammatic notation this 
simplifies the expansion of as developed for the 

Hubbard model in Ref. [33[ and, more recently, for the 
Bose-Hubbard model within an ITF in Ref. Q. With 
minor modifications, these diagrammatic rules can be 
taken over to the Schwinger-Keldysh formulation as fol- 
lows: 

1 . The building blocks of the expansion are the unper- 
turbed 2n-point cumulants. They are represented 
by a black circle located at site i with n in-going and 
n out-going legs. These legs carry the correspond- 
ing time variables. In-going legs are associated with 
a particle creation process, while out-going legs rep- 
resent a particle annihilation process. For example 
wc have for n = 1: 



Gj(ti; i 2 ). 



(24) 



2. A hopping process from site i to j is described by 
linking one out-going leg of the cumulant at site i 



6 



to an in-going leg of the cumulant sitting on the 
neighboring site j. We have to associate these in- 
ternal lines with a factor J. The number of internal 
lines within one diagram constitutes its hopping or- 
der. The time variable of an internal line has to be 
integrated. A subtlety coming from the CTPF in 
Keldysh space is the a 1 matrix between the vec- 
tor operators in Eq. (|17[) . which defines the matrix 
structure of the internal lines [3l| . 

3. Finally, the remaining legs have to be closed by a 
current, which are represented by a black square: 



reads: 



(25) 



The distinction between currents j and conjugate 
currents j is taken care of by the direction of the 
line, since j (j) is related to creation (annihilation) 
processes. Time and space variables of the currents 
have to agree with the ones of the corresponding leg 
and cumulant. 

In the free energy, all variables of any diagram have 
to be summed or integrated, so we can save space 
by completely suppressing the variables as well as 
the sums and integrals. For instance, the following 
closed graph has to be interpreted as 



E 



x s ij j i (ti)Ci(ti,t2)ij(t 2 )- 



dt 2 



(26) 



(27) 



As already mentioned, we will obtain a Ginzburg- 
Landau functional by performing a Legendre transforma- 
tion. We therefore note that the order parameter fields 
and the currents are conjugate variables: 



<a 4 (t)) Ho = *,-(<) 



(28) 



It might be astonishing that also the physical observables 



(29) 



which naturally cannot take different values at the same 
time, appear as contour-ordered quantities. However, 
we will later find equations of motion which determine 
the order parameter fields (t) and we will see that 
the ansatz tyf (t) = *&^(t) yields the physical solution 
of these equations. Nevertheless, for the time being, we 
have to keep the path index P during the envisioned Leg- 
endre transformation. 



5. Up to a given order in the hopping and in the cur- 
rents, the free energy is found by writing down 
all topologically inequivalent linked diagrams. In 
this diagrammatic notation, the statement of the 
linked-cluster theorem becomes literal: Diagrams 
of distinct cumulants not being linked via a hop- 
ping process do not contribute. The only remaining 
subtlety is the weight of each contributing diagram: 
On the one hand, we get a factor l/(mj!m 3 !mj!) 
in the mjth hopping order and 2m,, th order in the 
currents from the Taylor expansion of (|16[) . On 
the other hand, within one diagram there are ex- 
actly mjlrrijlmjl lines, which can be interchanged. 
If each permutation would lead to a different term 
within the cumulant decomposition, both factors 
would cancel. But as we will explicitly see in the 
fourth-order terms below, interchanging two lines 
or two vertices, might give exactly the same dia- 
gram again. To avoid an over-counting of these di- 
agrams we have to divide them by their respective 
symmetry factor, i.e. by the number M of possible 
permutations yielding the same diagram. 

Setting h = 1 from now on, the whole expansion up 
to first hopping order and fourth order in the current 



IV. RESUMMATION VIA LEGENDRE 
TRANSFORMATION 



We now define the Legendre transformation in the 
standard way, yielding a functional T[if, that we will 
refer to as the effective action of the system: 

/CO 
dt^-^+C.C.]. 
-oc 

(30) 



To obtain an expansion of T as a power series in J, \& and 
we first insert the expansion of J-[j, j] from Eq. (|2T|) 
into the definition of in Eq. (|28|) . yielding an expres- 
sion of the order parameter fields as a power series in J, 
j, and j. In the diagrammatic notation, the derivative in 
Eq. (|2"5|) is obtained by taking away a black square with 
an in-going leg from the graphs in T and, if there are 
more than one such squares within one graph, by apply- 
ing the usual product rule of differentiation. With this 
we obtain for the order parameter fields, which we denote 
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by white squares, the following expression: 



*i(t) = 



(31) 



Inserting the result from Eq. (|36l) into Eq. (|30j) together 
with (|27l) and discarding any diagram, which is higher 
than first order in the hopping or higher than fourth order 
in the order parameter fields, we finally find the effective 
action in the shape of a Ginzburg-Landau functional in 
first hopping order: 



-D*00 



(37) 



Variables fixed by the derivative appear explicitly in the 
diagrams, for all other variables the previous summation 
convention continues to hold. 

Inverting expression (1311) iteratively in J and j, we find 
a power series of \&](J). In zeroth hopping order 
and first order in the currents, only the first diagram of 
Eq. pip must be taken into account. Thus, the inversion 
reads: 



(32) 



where the white circle stands for the inverse cumulant 
[Cj(ti;t2)p . Note that the triangular structure of 
Ci(ti;t 2 ) is conserved under inversion: 



[C i (t 1 ;t 2 )} 



(33) 



Cf A{0) (t i; t 2 ) 



AE(O) 



EE(O) 



(*i;*2) 



where we have introduced the abbreviation 
-l 



C 



EE(O) 



(fl,*2) 



(34) 



-1 r 



c 



AE(O) 



(*1,*2) 



Re- inserting this in the second diagram of Eq. (13"TT) . we 
get the inversion up to the first hopping order: 



(35) 



and, finally, taking into account the currents up to third 
order, we get 




(36) 



Note that the lines between a cumulant and its inverse 
are no internal lines, i.e. they don't represent a hopping. 



/oo 
- oo 




Comparing the graphical content of the effective action 
in Eq. (|37|) with the corresponding one of the free energy 
in Eq. (|27[) . we find that it has a much simpler structure: 
In the last two diagrams of the free energy the hopping 
connects two cumulants, yielding a so-called one-particle 
reducible diagram, i.e. a diagram which might be di- 
vided into two diagrams by cutting one single line. In 
the effective action, however, these diagrams do not ap- 
pear and only the one-particle irreducible diagrams re- 
main [II IH. 

At this place, it is helpful to express the result (|37p in 
analytic terms. However, while in the diagrammatic ap- 
proach, apart from a straight-forward re-definition of the 
diagrams, it makes no difference if we work in real space 
or in Fourier space, the analytic expression is simplified 
considerably by a Fourier transformation. Therefore, we 
define the transformation into frequency space 



and into wave vector space 



dt f{t)e il 



(38) 



(39) 



For conjugate variables, e.g. *ff, the corresponding con- 
jugate transformations hold. 

Since the unperturbed Hamiltonian Ho in (|13[) is not 
explicitly time-dependent, the unperturbed cumulants 
may only depend on the differences in their time vari- 
ables. As we will explicitly see in the Appendix [A] where 
all the relevant cumulants are calculated in frequency 
space, this independence from an absolute time yields 
a Dirac <5-function in frequency space. Thus we are able 
to define 



(40) 



and we have 



C(a>i, uj 2 ; w 3 , w 4 ) — S(u>i + uj 2 - w 3 - cj 4 ) . (41) 

Regarding the spatial variables, the homogeneity of J,j 
assures that the hopping term becomes local by a trans- 
formation into wave vector space. In a cubic lattice with 
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lattice spacing a, for instance, we get 

3 



Jij — > 2J^cos(fc l a) = Jjr. 



(42) 



As all cumulants are local, the transformation into wave 



vector space always yields a fc-dependence in form of Sj: g. 

Since also a summation over k must be performed, we 
can completely suppress the wave vector indices of the 
cumulants. 

The effective action then reads in terms of Fourier 
transformed quantities: 



r (4,D = _£ j d0J ^{[cHi^-oijgj^j+ix; / d M % 1+fe - 2 - fe -3- fe - 4 

k {k} 

x (wi) [c(wi)] _1 (w 2 ) [c(w 2 )] _1 C(wi, lo 2 ; w 3 , w 4 ) [c(w 3 )] _1 *£ 3 M [cMP 1 *£ 4 (w 4 ). (43) 

I 



Now we can show that this expression of the effective 
action contains resummed hopping processes. To this 
end, we consider the first term in Eq. (|43|) . which is 
of second order in the fields: [c(oj)] _1 [l — J^a 1 c(uj)\ . 
Expanding this expression into a geometric sum, yields 

( c ( w ) S n [^k fjlc ( UJ )] n ) ■ Diagrammatically, this rep- 
resents the inverse of a sum of arbitrarily long hopping 
"chains" . Thus, by means of a resummation, any of these 
terms is automatically taken into account within the ef- 
fective action. This is the clue which transforms the ini- 
tial perturbative hopping expansion for the free energy 
into an effective 1/D expansion for the effective action, 
from which we may expect reliable results not only in the 
MI but also in the SF phase. 

We conclude this section by showing how to extract 
physical information from the effective action given by 
Eq. (14*3")) . The easiest and most direct way to do so is by 
noting that the physical situation corresponds to vanish- 
ing currents. Due to Legendre identities this means 



ST 



0. 



(44) 



refer to these equations as the equations of motion as 
they determine the order parameter field. In the next 
section we will deal with their solutions. This will yield 
both the phase boundary, if we consider the static case, 
and the excitation spectra in the dynamical case. 



V. EQUATIONS OF MOTION 

Analyzing the structure of the effective action T in 
(|43f we find that the second-order term ^ E \I/ E is re- 
lated to the vanishing cumulant C AA and the fourth- 
order term ij/ s \I> s \I> E \I> E is related to the vanishing cu- 
mulant C AAAA . This means that all non- vanishing terms 
in Eq. (|43|) contain ^ A at least to first order. Since the 
first component of Eq. (|4"4")) is given by ST/S^ = 0, all 



and the complex conjugate of this equation. We will 



the terms in this equation will contain 1 J fA or W as a 
factor. As we aim at solving the equations of motion with 
the ansatz 1 & + (t) — *&~(t), i.e. with an order parameter 
field which is independent of the time path, we are look- 
ing for solutions \& A = 0. Thus, we find that the first 
component of Eq. (|4"4")) is always trivially fulfilled. The 
second component reads: 



J 



ST 



pOO pOO pOO 

- dujt du 2 / du 3 k 1+ k 2 

^ J — OO J DC J — OO - r* r- 



-k^—k 



k± ,k 2 



C AAAS (°)( Wl ,u, 2 ;a, 3) u;) + C AASA (°)( Wl , W2 ; w ,w 3 )' 



*f^o*f>)*f> 3 ) 



G A (°) (wi)G A (°) (w 2 )G R (°) (w 3 )G A (°) (w) 



= 0. 



(45) 



The four-point cumulants C KiK2K3K4(0) (cji,u; 2 ;cj3,cj 4 ) 

are symmetric in {wi,Ki} -h- {cj 2 ,K 2 } and {w 3 ,K3} •<->• 

{w 4 ,K 4 }, such that C AAAE (°>(wi, w 2 ; w 3 , w 4 ) = 
c .aaea(o) ( Wi ( W2 . Thus ^ from the initially i 6 dif- 



ferent 4-point cumulants, only one has to be calculated. 
This is done in the appendix, where also explicit ex- 
pressions for the retarded and advanced Green functions 
qR(o) an d Q A (o) are derived. There we will see fur- 
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thermore that 2G AASA (°)(wi, w 2 , oj 3 , cj 4 ) coincides with 
the advanced four-point function C A ^ (uji,uj2', 0)3,014) 
which is defined as the thermal average of multiple 
commutators times Heaviside step functions (cf. |23|). 
As the complex conjugate of advanced Green functions 
are retarded Green functions, the complex conjugate 
of Eq. (|44l) amounts to interchanging the respective 
retarded and advanced functions. 

Due to the triple frequency integrals and the triple 
wave vector sums, Eq. (|4"5)l is not straightforwardly 
solved. In order to continue analytically, we now assume 
the equilibrium ansatz = S(ui)5j: of a spa- 

tially homogeneous and temporally constant order pa- 
rameter field. With that ansatz, all integrals and sums 



become trivial, and we get the algebraic equation: 

1 G A ( >(0,0;0,0)^ c s q 

G A (°)(0) 3 4G A ( )(0)G A (°)(0)G R ( )(0)G A (°)(0) ' 

(46) 

Later on we will read off from this equation physically 
important results like the phase boundary. At the mo- 
ment we only need the information that this equation 
defines an equilibrium value around which we per- 
form a harmonic approximation for the effective action 
Eq. (gSJ) with the ansatz *g (w) = + $f (w). This 
yields 



dwi 



2 , 

ki,k 2 

' ,5*5 (wi)(5* A (wa) 



duio 



S 2 T 



SVf (wa) 



S^(W 2 )!P|(W 1 ) + 



oq 

5 2 r 



s 2 r 



^|K)^ A ( W2 ) 



(47) 



cq 



cq 



<J*5 (wi)«*f (wa) 



cq 



Other derivatives either vanish or lead to terms which are 
of second order in \P A . They would vanish when the re- 
sulting equations of motion are solved by <J/ A = 0. Eval- 
uating explicitly the respective derivatives in Eq. (1471) . 
yields the linearized expression: 



G R (°)(w) 



Jt 



l*c 



: G R (°)(w,0;0,w) 



2G A (°) (0)G R (°) (0)G R (°) (w)G R (w) 



l*c 



; G R (°)( W ,- W ;0,0) 



fc v-/ 4 G a(o)(_ w ) G r(o)( ) G R(o)( w ) Gr(o) ( ) 

x^N = 0. (48) 

Now we will investigate how this linearized equation of 
motion is related to the Green functions. To this end we 
start with the statement 



Making use of the product rule for functional derivatives, 
we find the following identity 



5]t 



k 



A 



(50) 



and a similar one from = Sj A / <5j 
identity 



A . With the Legendre 



<5 J" 



(51) 



5 2 T 



5j?5jf 



_ G A 



J=j=0 



which is dual to (|4"4")l , we are finally able to relate second 

(49) 

derivatives of the effective action T to the advanced two- 
point functions: 



<5* A £ (-w)(5*^ g (-w) 



cq 



(5 2 r 




a 2 r 






<5* A g (-w)5* E £ (-w) 


eql 5*£(u;),5*£( W ) 




eq 5*f( w )5^ s (- w ) 


eq 



(52) 



and the complex conjugate expression for the retarded solved, when the superfluid Green function (|52p diverges. 
Green function. Note that the equation of motion (|48|) is 
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For J < Jpb at a given fi and U, only \& e 



solves 









2 
1 




-3 -2 - 1 

-1 


1 2 3 







FIG. 2: Spectral weights at the lobe tip (up) and for J/U 
0.02 (below) as a function of k = k • e x . 



However, the Green functions contain more information 
than the equation of motion, since they allow, in princi- 
ple, to extract also the respective spectral weights from 
its imaginary part as we will see below. 



VI. RESULTS 

The non-trivial solutions of Eq. (|4"5j) can be divided 
into two classes which we work out at T = 0. How to 
perform the T = limit for the respective Green func- 
tions and cumulants is shown in the appendix. 



Eq. (|46p . and thus we are in the disordered, i.e. insulat- 
ing phase. For J > Jpb, the system is superfluid. A plot 
of the phase boundary ([55)) for T — and < fi/U < 1 
can be found on the left side of Fig. [3] for n = 1. In 
the considered first hopping order, this quantum phase 
diagram is identical to the mean-field result Q with a de- 
viation from recent high-precision Monte Carlo data [36[ 
of about 25 %. Within a Landau expansion the second- 
order hopping contribution has recently been calculated 
analytically in Ref. [13] decreasing the error to less than 
2 %. A numerical evaluation of higher hopping orders 
has even been shown to converge to a quantum phase di- 
agram which is indistinguishable from the Monte-Carlo 
result [l3-[l7|. 

For obtaining dynamic solutions we have to consider 
the linearized equation of motion Eq. (1481) . Numerically, 
we obtain the dispersion relations shown in Fig. |3l where 
the MI-SF transition can clearly be observed by a quali- 
tative change in the excitations of the system. 



B. MI spectra 

In the MI phase, we find the gapped particle/hole exci- 
tations from mean-field theory [9[. At the phase bound- 
ary, we have to distinguish the tip from the rest of the 
lobe in accordance with the critical theory of the BH 
model [111 [3O] - While at the tip of the mth lobe, i.e. 
at [imjU = y m(m + 1) — 1, both excitations become 
gapless and linear for small |fe|, off the tip at /J, > /i m 
(fi < fj, m ) only the particle (hole) mode becomes gapless 
and remains with a finite effective mass. 

Since in the MI phase and at the phase boundary we 
have I* 



cq 



0, we are able to determine analytic solu- 



tions of the equations of motion Eq. (|4"5|) : 



1 

2h 



2/i 



■k 



U(2n- 

± J IP- 2JrU(2n+l) 



(54) 



The Green function in Eq. (|52|) reduces in the MI 
phase to 



A. Phase boundary 

Static solutions are obtained from Eq. (|4"oT) . The right- 
hand side of Eq. (|46|) does not explicitly depend on the 
hopping parameter, but it is multiplied by I'J'oql 2 > 0. 
Thus, its sign is independent from J. The left-hand side 
of Eq. (|4"o) depends explicitly of J and changes its sign 
at some critical value JpB(fJ-,U), which is obtained by 
setting the left-hand side of Eq. (J46J) equal to zero. From 
(|A3]) and (|A9|) we get explicitly 



Jpb(m, U) = - 



(Un- fj,-U)(Un- /x) 
6(/VC/+l) ' 



(53) 



G R(MI) 



1 



GK(°)(w) J ^ 



(55) 



Comparing this with the general spectral representation 
of the retarded/advanced Green functions [33] 



du/p(o/, k) 



— ± iirS^LO — lo') 



(56) 



where V denotes the principal value of the integration 
across the singularity and p(u, k) denotes the spectral 
function, we can also obtain information about the spec- 
tral weights of the excitations, i.e. how much of the total 
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FIG. 3: Excitations in k — (1,0,0) direction are plotted in Al - CI for different values fi/U and J/U, which are marked in 
the phase diagram (left). In the MI phase (green) and on the phase boundary (blue), the two T — modes can be interpreted 
as particle (dotted) and hole (dashed) excitations. At the tip of the Mott lobe (B), both modes become gapless, whereas for 
larger (smaller) chemical potentials fj, only the gap of the particle (hole) mode vanishes. In the SF phase (red), the gapless 
mode turns into a sound mode, but also a gapped mode is present everywhere in the SF phase. The smooth transition from 
the MI excitation to the SF excitation is further analyzed in A 2 - C 2 and A 3 - C 3, where the effective mass m and the 
gap A of each mode are plotted as a function of J/U. The sound velocity c of the massless SF excitation, plotted in A3 - C3, 
vanishes at the phase boundary except at the tip, indicating the existence of a different universality class in this configuration. 



excitation energy is stored in each excitation. To this end 
we have to divide Eq. (1551) into its real and imaginary 
part, where the latter is due to an infinitesimal ie-shift of 
the poles from the real axis into the complex plane (see 
Appendix A). By performing the limit e — ¥ 0, we find 

p(ui, k) = 5{uj - uj + (k))w + (k) + S(uj - u-(k))w-(jt) (57) 
with the respective weights 

,a 1 / U(l + 2n)-Jr \ 

w ± (k) = -\l± K I k (58) 

2 \ ^-2J s U(2n + l) + Jl) 

Note that they do not depend on the chemical potential. 
At the tip of the lobe, both weights diverge at k — 0. We 
can check the spectral function ([ST|) . (|58p by noting that 



it obeys the sum rule [37j: 

/oo 
dwp(u,k) = l. (59) 
-oo 

The spectral weights are plotted for two different ratios 
J/U in Fig. H 



C. SF spectra 

Turning into the SF phase, the gapless mode rapidly 
looses its mass and has to be identified with the Gold- 
stone mode which arises due to the broken U(l) symme- 
try. Indeed, within the Ginzburg-Landau theory it turns 
out in the limit k — > and u) — > 0, that the excitation 
(5^/(oj, fc) stems from variations of the phase. Within a 
slave-boson approach it has even been shown in Ref. [lfjj 
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that, also for general wave vectors k, phase variations 
dominate this excitation. Thus density variations arise 
which make this mode sensitive to Bragg spectroscopy. 
Recently, the whole sound mode has been measured via 
Bragg spectroscopy far away from the phase boundary 
and could be well described via a Bogoliubov fit @ ■ An- 
alytical results for this mode can be obtained by expand- 
ing Eq. (|4"8")l for small uj and small k up to second order. 
The resulting algebraic equation is solved by the ansatz 
u> = c\k\. Defining the dimensionless quantities jl = /i/U 
and J = J/U, we find that the sound velocity c of the 



Goldstone mode at T 
sionless quantity 



0, which is given by the dimen- 



c(n,fl, J) 



z(n,iJ,,J,U) 
all/h 



(60) 



along any lattice direction. Here n gives the occupation 
number of the MI lobe, above which the SF theory is 
constructed. For n = 1, for instance, we find the explicit 
result 



c(l, ft, J) = | [j/i 2 (l + /i) 3 (3 - 2/2)( - 3 + 8£ - 10fl 2 + 4fi 3 + jl 4 ) 2 {(fl - l)/x + 6J(1 + /}))] / 



36J(1 - /}) a /} a (-27 + 108/i + 9/2 2 - 92/2 3 + 3/i 4 - 24/i 5 + 7/t B ) - (-1 + /i) 3 /2 3 



(27 - 135/i + 36/i 2 + 172^ - 210// 4 + 294/t 5 - 196/i 6 + 60£ 7 
(27 - 270/i + 1359/i 2 - 3860/i 3 + 5950/i 4 - 4512/2 5 + 1198/i 6 + 100/i 7 

I 



15£ 8 + m 9 ) 



135/t 



(61) 



18 J 2 (1 + /i) 2 

5 - 66/i 9 + 3^ 10 )] } 



1/2 



Additionally to that sound mode, however, also a gapped 
mode survives the quantum phase transition or arises 
again if we depart from the lobe tip. It can be smoothly 
mapped onto one of the respective MI modes, which is 
shown in the plot of both the effective masses and the 
gaps on the right side of Fig. [3J The existence of such a SF 
gapped mode is in accordance with theoretical results ob- 
tained previously in Refs. [1CH12| .and it has also recently 
been confirmed experimentally [7j. This mode is inter- 
preted in Refs. [3, E3] as an amplitude excitation which 
corresponds to an exchange between condensed and non- 
condensed particles at constant overall density. Although 
Eq. (|48[) does not allow pure amplitude excitations, we 
can back this interpretation of predominant amplitude 
excitations by observing that the zero-momentum energy 
transfer at the phase boundary corresponds to the cre- 
ation of a particle/hole pair. 



D. Critical exponents 

From our excitation spectra we can directly read off 
the dynamical critical exponent z, which is z = 1 at the 
tip, where a relativistic dispersion relation is found, and 
z = 2 off the tip, where the system behaves like a normal 
Bose gas. This is in agreement with the critical theory 
of the BH model [2!| [3(| ■ Physically, the difference be- 
tween the "generic" transition off the lobe tip and the 
"XY-like" one at the tip is caused by the density varia- 
tions that are only absent at the lobe tip. Writing the gap 
A as a function of J — Jpb(^), i.e. of the distance to the 
phase boundary, we find that A ~ [J— Jpb(h)] zl ' with the 
mean-field exponent v = 1/2 for any vertical phase tran- 
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transition : 
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/ /, / XY transition : n 






J z = 1 . V = 1/2 
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/ // generic transitions: 






'// z = 2 , V = 1/2 
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■ 





0.5 1 

FIG. 4: The critical exponents depend on position and di- 
rection of the phase transition. The grey lines of constant 
densities in the SF phase are obtained via a thermodynamic 
effective action based on an analogue expansion in imaginary 
time 0. 



sition, see Fig. 2) Since vertical transitions really cross 
the phase boundary, their exponents can be obtained an- 
alytically by considering only the behavior within the 
Mott phase. At the tip, there is also the possibility to 
touch the phase boundary horizontally, i.e. without en- 
tering the MI phase. We find that this is a special case 
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of the generic transitions with zv = 1, where z = 1 and 
^ = 1 are obtained numerically. 



E. Gross-Pitaevskii limit 

Surprisingly, also the regime deep in the SF phase turns 
out to be accessible with the Ginzburg-Landau theory. 
Expanding the Green functions in U, in the lowest non- 
trivial order they do not depend on temperature and re- 
duce to 

G R/A(0) (wi,w 2 ) = -2tt— !— <S(wi - w 2 ) + 0(U) (62) 
/-t + oj 

C R / A W(wi,W2; lo 3 , Ui) = —2it5(ijJi + lo 2 — w 3 — W4) 
2U 



0(U 2 



(63) 



(^ + uji)((i + w 2 )(/i + wa)(p + w 4 ) 

Inserting this into Eq. (|45|) leads, after a Fourier trans- 
formation into real time and space, to the equation of 
motion 



which is the lattice version of the Gross-Pitaevskii (GP) 
equation (38|. From this follows the Bogoliubov sound 
mode of a fully condensed system [9( , which is given by 



w fc = \/ e k + 2Une k> 



(65) 



where = 6 J— denotes the free dispersion. For small 
\k\, this dispersion is linear with the sound velocity 



—- = y/2J((l + 6J), 

a/n 



(66) 



which coincides with the limit ell from Eq. (|61[) for 
U -> 0. 

Since in this limit the sound mode predicted by our the- 
ory agrees perfectly with the one derived from the GP 
equation, it is interesting to inspect also the gapped mode 
which is not predicted by the GP theory. Thus, we can- 
not first expand the equation of motion for small U and 
then consider the solutions, since this would reduce our 
equation of motion to the GP equation. Instead, we 
may investigate what happens to the solutions of our full 
equations of motion for small U. Numerically we find 
that the gapped mode converges to the constant disper- 
sion u(k) = 2/i, if U becomes negligibly small. As the 
chemical potential \i corresponds to the energy needed 
to add or take away one non-interacting particle from 
a lattice site, this excitation may be interpreted as the 
creation of a particle-hole pair. 



VII. COMPARISON WITH EQUILIBRIUM 
THEORY 



Finally we compare our results with the ones obtained 
within a similar Ginzburg-Landau theory in imaginary 



time [14| . Whereas at T = both formalisms yield iden- 
tical results, a mismatch occurs for finite temperature. 
This is surprising, since both formalisms are considered 
to be equivalent in equilibrium (23l . l2~ij , but there also ex- 
ists a number of papers, in which possible di sag reements 
of the CTPF and the ITF are discussed [H-E3. 

In order to localize this mismatch, we note that the 
equations of motions in the ITF have the same structure 
as ours in Eq. (|48|) . but the retarded and advanced cumu- 
lants have to be replaced by the analytical continuation of 
the thermal Green functions to real frequencies. For the 
2-point function we find that the analytical continuation 
yields exactly the retarded/advanced functions. How- 
ever, in case of the 4-point functions, the thermal Green 
function has a decomposition into lower cumulants, while 
the retarded and advanced Green functions do not have 
a similar decomposition (see Appendix A and Ref. [HI). 
Apart from this missing cumulant decomposition terms, 
both formalisms agree completely, but these terms, which 
vanish for T = 0, might become large if temperature in- 
creases. 

We note that the decomposition vanishes for retarded 
and advanced Green functions due to the choice of our 
time-evolution contour in Section |TTJ On the Keldysh 
contour the path of the largest-time operator plays no 
role, which results in the vanishing of contour-ordered 
products of A-indexed operators and therefore of any 
possible decompositions of the retarded and advanced 
Green functions. Obviously, the path of the largest-time 
operator would play a role, if not both forward and back- 
ward path were chosen to be along the real time axis, e.g. 
if we shifted the backward axis by — i/3 as proposed in 
Ref. [44( . Thus we can conclude that the Keldysh formal- 
ism working with a purely real time-evolution contour is 
not able to produce the correct equilibrium configuration 
of the full system for finite temperature. This result is 
supported by Ref. [45| , where it is argued that throwing 
away the imaginary part of the time-evolution contour 
reduces the generating functional Z to 1, whereas along 
the full contour it really represents the partition function. 

These formal argumentations can be backed by 
the physical picture that an instantaneous change of 
one time-independent Hamiltonian to another time- 
independent Hamiltonian should lead to the relaxation 
into a new equilibrium with a new temperature. If we 
neglect the imaginary part of the time-evolution contour, 
we formally replace the thermal average with respect to 
the new Hamiltonian by a thermal average with respect 
to the old, unperturbed Hamiltonian, without modifying 
the temperature. This might be a justified approximation 
in a non-equilibrium system where temperature can only 
be defined in the unperturbed initial state, in our time- 
independent system, however, it yields wrong results. 

The agreement of both formalisms at T = can be un- 
derstood as a consequence of the Gell-Mann-Low theorem 
22] stating that the systems remains in the ground-state, 
if a perturbation is adiabatically switched on. Since the 
Keldysh ansatz has pushed this switching into the infi- 
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nite past, no additional assumptions about its adiabatic 
properties had to be made (24[. 



VIII. SUMMARY AND OUTLOOK 

We have developed a real-time and finite-temperature 
Ginzburg-Landau theory for the Bosc-Hubbard model by 
applying the Schwinger-Keldysh formalism of a closed 
time-path reviewed in Section HU The perturbative hop- 
ping contributions for the free energy have been re- 
summed via a Legendre transformation, yielding a large 
D expansion for the effective action. Therefore, we have 
been able not only to extract the phase boundary from 
the equations of motion (|4*4l . but also to calculate the 
excitation spectra in the SF phase: There, the parti- 
cle/hole excitations from the MI phase turn smoothly 
into a gapped amplitude and a gapless phase mode. Sur- 
prisingly, we even obtain reliable results in the limit of 
small U, where our theory turns into the Gross-Pitaevskii 
theory. 

A comparison with a similar theory making use of 
imaginary times (l4j . has shown an unexpected mis- 
match due to the negligence of the imaginary part of 
the time-evolution contour. In the future they should 
be overcome by a new calculation, which explicitly in- 
cludes this imaginary path. At T = 0, our theory al- 
ready works in full agreement with ITF and could, thus, 
be directly applied to Bose-Hubbard Hamiltonians with 
time-dependent hopping parameters J, which are rele- 
vant in collapse and revival experiments [l8]. Including 
a possible time-dependence of U might be more prob- 
lematic, but since the relevant quantity is J/U, we are 
be able to keep U constant and put effectively all time- 
dependencies into the hopping parameter [4f| . Then the 
only difficulty consists in solving the equation of motion 
which then cannot be linearized around an equilibrium 
any longer. By applying numerical methods, one should 
be able to obtain non-equilibrium results with the for- 
malism developed here. 
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Appendix A: Calculation of the Cumulants 

As worked out in Section [Vj the equations of motion 
depend on the retarded and advanced 2- and 4-point 
cumulants. Written in the Keldysh basis, the retarded 
and advanced Green functions are given as averages of 
contour-ordered operator products with one operator s 



and the rest being operators 6 A . If the S-indexed op- 
erator is an annihilation operator, we have the retarded 
function, otherwise it is the advanced one. From this, 
it can be directly seen that advanced and retarded func- 
tions are linked via complex conjugation. We thus need 
to calculate explicitly only one of both. Any cumu- 
lant decomposition of these functions necessarily involves 
(T c O A A ) or (0 A ), both being zero, since they contain 
exclusively A-indexed operators. Therefore, the decom- 
position of the retarded/advanced Green functions con- 
tains Green functions of lower order, and we don't have 
to distinguish between retarded/advanced cumulants and 
the corresponding Green functions. However we should 
note that the unperturbed cumulants are always local, 
whereas the unperturbed Green functions may also de- 
scribe independent processes on distinct sites. In this 
appendix, however, we completely suppress spatial vari- 
ables and assume locality for all objects. 

Starting with the 2-point function, we first re-write it 
by making use of the Heaviside step function: 

G R (°)(ii;i 2 ) =z(f c a s (i 1 )a A t(t 2 )\^ (Al) 

= ie{t x - t 2 )Tr[e-' sSo a(t 1 )at(t 2 ) - ^(h)^)] /2 (0) , 

where — Tr e _ ^ H °. The traces are best calculated 
in the occupation number basis solving the eigenvalue 
problem 

H |n) = E n \n) , (A2) 
where the energy eigenvalues are given by 

E n = —n(n - 1) - pm. (A3) 
The retarded Green functions can thus be written as 

oo 

G R(0) (h,t 2 ) = i9(h ~ t 2 ) e-^ \(n + 1) (A4) 

/Z (0) 

with the partition function 



xe -i(E„ +1 -E„)(t 1 -t 2 )_ ne 



n=0 

j(S„-E„_ 1 )(t 1 -t 2 ) 



0E„ 



(A5) 



n=Q 



In order to determine the Fourier transform of (|A4I) ac- 
cording to Eq. (|38p , we use the integral representation of 
the step function: 



t 2 ) = lim 



x + ie 



(A6) 



In the following we will suppress the limit-symbol for 
reasons of brevity. Furthermore, we will make use of the 
Fourier representation of the Dirac (5-function: 



6(lui — LU 2 ) 



I 

2^ 



(A7) 
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Performing the substitution t = ti — we get at first 



n=0 



Ax 

-oo^ + ie /_ 



at 



-ne 



-i(E n -E n -i-u 2 +hx)t 



(A8) 



Evaluating both integrals yields then the result 

1 oo 

n+ 1 



n=0 
7? 



■En+i — -En — ^2 — ie E n ~ E n -\ — u>2 — ie 



.(A9) 



To get the advanced 2-point function, we must only re- 
place the minus sign in front of the e by a plus sign. 

The procedure for calculating the retarded 4-point 
function is much the same. At first we express it in terms 
of Heaviside functions: 



C 



EAAA 



6(ti - h)9(t 2 - t 3 )9(t 3 - U) 



a(ii),a(t 2 ) ,a f (i 3 ) 



x6(t3-t 2 )0{h-t 4 )( 



a f (i 4 ) 



+ 0(ti-t a ) 



a(ti),a t (i 3 ) ,a(i 2 ) 



-6(h - ta)0(t 3 - t 4 )0(t 4 - 1 2 ) 



,a(ta) 



H ' t 3 +±t4 



(A10) 
I 



Here, the symbol t% O (4 means that we still have to 
symmetrize the expression in these variables. Apart from 
a factor 2, it is identical to the usual definition of retarded 
n-point functions [39| . We define 



2C 



EAAA 



(il,^2;^37^4) = C R (^1, £2! ^3 ; ^4) 



(All) 



Again the thermal averages can be evaluated by tracing 
the operator products in the occupation number basis. 
Then the resulting expression is Fourier transformed in 
the same way as before when dealing with the 2-point 
function. Since the number of different terms is much big- 
ger now [43[ , we do not write them down explicitly. The 
whole function depends only on time differences, thus 
one Fourier transformation leads to the Dirac function 
5{u)i + u>2 — 0J3 — uji). From the Fourier representation of 
the three Heaviside functions, we now get three different 
infinitesimal frequency shifts iiej. The real part of the 4- 
point function is found by setting all ei = 0, whereas the 
imaginary part is obtained by carefully taking the limits 
€i — > 0. Since the knowledge of the latter is only needed 
for the spectral weights but not for the spectrum itself, 
we only give here the real part of the 4-point function, 
which has a relatively compact form: 



Re C R (wi, L02; u>3, W4) = — 2tt5(lui + L02 — 0J3 — UJ4) /J 
j n(n+l) 



-0E„ 



n=0 
1 



z(°) 



{E n — E n -l — Uli)(E n +l — E n — LO4) \E n — E n -1 — LO3 E n +l — E n — LO2 

n(n+l) ( 1 1 



+ 



+ 



(E n — E n -1 — 0J3)(E n +l — E n ~ LO\) \E n — E n -1 — L0 2 E n +l — E n — LO4 

-(n + 2)(n + l) 



(E n +1 — E n — L04)(E n +2 — E n — U>\ — LO2) \E n +l — E n — LOl E n +1 — E n — LO2 

n(n — 1) 



(E n — E n -1 — ljJ4)(E n — E n -2 — W\ — ^2) \ E n — E n -\ — U3\ E n — E n -\ — LO2 

2(n+l) 2 

(E„+i — E„ — LOi)(E n +i — E n — u>2)(E n +i — E n — W4) 
-2n 2 ) 



(E n — E n -i — u)i)(E n — E n -i — bj 2 )(E n — E n -\ — L04) 



(A12) 



1G 



If we wish to take the zero-temperature limit j3 — > oo 
in Eq. (IA9|) or (|A12|1 . we must note that there is one 
occupation number no, which represents the commensu- 
rate ground-state of the unperturbed system, and thus 
E no < E n for any integer n ^ hq. By factoring out 



e -/3E„ - m ^ e B itzmann sums both in the denominator 
and the numerator, all terms in the sum with n =/= hq 
remain with a factor Q-Pi E ^ E n a ) q f or p O0i Thus 
only the term with n = hq is not suppressed and survives 
the zero-temperature limit. 
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